Stellar configurations in f{R) theories of gravity 
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We study stellar configurations and the space-time around them in metric f{R) theories of gravity. 
In particular, we focus on the polytropic model of the Sun in the f{R) = R~ /R model. We show 
how the stellar configuration in the f{R) theory can, by appropriate initial conditions, be selected 
to be equal to that described by the Lane-Emden -equation and how a simple scaling relation 
exists between the solutions. We also derive the correct solution analytically near the center of 
the star in f[R) theory. Previous analytical and numerical results are confirmed, indicating that 
the space-time around the Sun is incompatible with Solar System constraints on the properties of 
gravity. Numerical work shows that stellar configurations, with a regular metric at the center, lead to 
fppN — 1/2 outside the star i.e. the Schwarzschild-de Sitter -space-time is not the correct vacuum 
solution for such configurations. Conversely, by selecting the Schwarzschild-de Sitter -metric as the 
outside solution, we find that the stellar configuration is unchanged but the metric is irregular at the 
center. The possibility of constructing a f{R) theory compatible with the Solar System experiments 
and possible new constraints arising from the radius-mass -relation of stellar objects is discussed. 
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I. INTRODUCTION 

Current cosmological observations provide strong ev- 
idence against a critical density matter dominated uni- 
verse. Observations on supernovae type la [l|, cosmic 
microwave background [2] and large scale structure [3| 
all indicate that the expansion of the universe is not pro- 
ceeding as predicted by general relativity (GR), if the uni- 
verse is homogeneous, spatially flat and filled with non- 
relativistic matter. The underlying assumptions have 
hence been questioned, resulting in the emergence of a 
cosmological concordance model, the ACDM -model, de- 
scribing a flat and homogeneous universe dominated by 
cold dark matter and dark energy in the form of a cosmo- 
logical constant. Alternatives to the simple cosmological 
constant are, however, numerous (for a review see e.q. 

and other assumptions, such as e.g. homogeneity [5| 
have been questioned. 

An alternative route to solving the dark energy prob- 
lem is to consider modifying the underlying theory of 
gravity i. e. relaxing the assumption that general relativ- 
ity is the correct theory of gravity on cosmological scales. 
A popular choice is the class of f{R) gravity models that 
has received much attention in the recent literature (see 
e.g. and references therein). 

In a /(i?) gravity model, deviations from general rel- 
ativity arise by generalizing the Einstcin-Hilbert -action 
with an arbitrary function of the curvature scalar, f{R). 
Such modification has to face many challenges that gen- 
eral relativity passes, including instabiliti es fill [l^. fivj . 
solar system constraints(see e.g. [l^ Hi l2Cll| and refer- 
ences therein) and evolution of large scale perturbations 
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In particular, the solar system observations offer a 
good testing ground for any modification of general rel- 
ativity by comparing the Parameterized Post-Newtonian 
(PPN) parameters [2^ |M, |25|, III] with observations. 
This question has recently been extensively reviewed and 
discussed by a number of authors [13, [H, [H, [s^l in order 
to determine the relevance of the Schwarzschild-de Sitter 
(SdS)-solution in the solar system (recently a new class 
of models that can evade the Solar System constraints 
has been introduced [3l[). The SdS-metric is an exact 
vacuum solution in a large class of /(i?)-theories of grav- 
ity that is in agreement with all solar system observa- 
tions with an appropriate cosmological limit. However, 
the higher order nature of f{R) theories makes the SdS- 
solution (see eg. [3^, |33|) non-unique. This property of 
f{R) theories also demonstrates itself on a cosmological 
scale, making reconstruction of the form of f{R) from 
the expansion history of the universe non- unique [34| . 

As a result of the recent discussion, it has become clear 
that the SdS metric is unlikely to be the correct vacuum 
solution in the Solar System for the 1/R model. Instead, 
the PPN Solar System constraints are valid in a limit 
that corresponds to the limit of light effective scalar in 
the equivalent scalar-tensor theory. This is equivalent to 
requiring that one can approximate the trace of the field 
equations by Laplace's equation [HI in the corresponding 
f{R) theory. This result has now also been considered 
by numerical calculations [ssj , where the field equations 
are integrated numerically from the center of a star for a 
fixed matter distribution. 

Relatedly, in a recent work [3^ we considered perfect 
static fiuid sphere solutions in f{R) theories of gravity. 
Again, the higher order nature of the f{R) gravity the- 
ories demonstrates itself in that the mass distribution 
alone does not uniquely determine the gravitational the- 
ory, unless the boundary conditions are fixed. If one im- 
poses the SdS-metric as a boundary condition, one finds 
that the solutions are constrained. 

Here we consider these questions by numerical and an- 
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alytical means. We solve the set of field equations both 
inwards and outwards, i.e. by starting from the center 
and the boundary of a star. In contrast to ^35|, we do 
not fix the mass distribution beforehand and for com- 
pleteness also study configurations with non-negligible 
pressure. Furthermore, we consider mass distributions 
with the SdS-metric as a boundary condition and discuss 
appropriate analytical limiting solutions near the origin. 



II. f[R) GRAVITY FORMALISM 

The action for f{R) gravity is 

j d'x^g(^f{R)+C„) (1) 

and the corresponding field equations derived by variat- 
ing wrt the metric g^y are 

F{R)R,,, - i/(i?).g^. - (V,,V, - g^.n)F(i?) = 8^GT^,. 

(2) 

Here T^i, is the standard minimally coupled stress-energy 
tensor and F{R) = df/dR. Contracting the field equa- 
tions we get another useful form: 

F{R)R - 2f{R) + 30F{R) = 87rG(p - 3p). (3) 

We consider spherically symmetric, static configurations 
{p = Pif)j P = pi''')) ^'^d adopt a metric: 



ds^ 



B{r)dt^ - A{r)dr'^ ~ r^{de^ + air? dcf^). (4) 



In this metric, the non-trivial component of the continu- 
ity equation D^T'^'^ — reads: 



P 



2B 



(P + P), 



(5) 



where ' = d/dr. Note that, like in GR, the equation 
of continuity is automatically satisfied [s^l and hence no 
additional information is gained on top of the field equa- 
tions. On the other hand, one can choose the continuity 
equation as one of the equations to be solved instead of 
using the full set of field equations. 



A. Equations 

For a given action function f{R), one can in principle 
take any suitable set of modified field equations, Eq. ([2]) 
along with the equation of state, p = p{p), and solve 
for p. A, B. However, this can in practice prove to be 



problematic, since the modified Einstein's equations are 
highly non-linear and high order differential equations, 
up to fourth order in B and third order in A. 

To simplify the problem, it is useful to consider f{R) 
and F{R) as independent functions of r. In order to cor- 
rectly account for their dependence, one then needs to 
supplement the set of field equations with an appropri- 
ate additional constraint, / — f{F), determined by the 
details of the particular f{R) theory in question. In this 
description we are able to view F{r), A{r), B{r) and p as 
the fundamental set of unknown functions to be solved. 
Note that unlike before, now the equation of continuity 
is not automatically satisfied but is an additional, inde- 
pendent differential equation. This is due to the fact that 
F(r) (and /(r)) is not given in terms of the scalar curva- 
ture but is viewed just as a independent function of the 
radial coordinate r. 

Thus the modified Einstein equation along with the 
equation of continuity forms the set of independent equa- 
tions to be solved. These nonlinear equations are, how- 
ever, only second order in F and B, and first order in 
A and p requiring in total six initial conditions for com- 
pletely determining the solution. This is apparently one 
less that is needed if one proceeds by solving the field 
equations directly in terms of A and B, demonstrating 
that the higher derivatives of the modified field equations 
only appear in the combinations of the derivatives of the 
scalar curvature R. 

In this paper, we consider the CDTT-modcl introduced 
in [6] with 
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but generalization to more complicated models is 
straightforward. In this particular case, it is easy to see 
that the equation relating / and F is 



f = p' 



F 



y/F^ 



(7) 



In a more general case, like f{R) = i? — p.'^/R. + R'^/f3'^, 
the functional relation is more involved and in general 
one may need to resort to numerical means. 

As is well known, the CDTT-model has a homogeneous 
solution of de Sitter -metric with constant scalar curva- 
ture R = Ro = —V^p"^. In order to have the desired 
physically plausible late time behaviour, we therefore set 

Using Eq. ([7]) and defining F{r) = 4/3 -I- v{r), n{r) = 
B'{r)/B{r), straightforward algebraic manipulations re- 
sult in the following set of equations: 
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This set of equations is suitable for numerical integration 
made in this article. 

Note, that the dependence on B{r) dissapears com- 
pletely from the equations ([5]), i.e. the equations are 
only first order in n(r) reflecting the leftover free time 
scaling of the metric component B(r). 



B. Solution near the origin 

In order to solve the field equations numerically from 
the center, one also has to consider the question of bound- 
ary conditions. Clearly, one cannot start the numerical 
integration from the origin r = due to singularities, but 
a small distance from it. Therefore one must select the 
initial values such that they correspond to a desired and a 
possible solution: otherwise, one might start the calcula- 
tion from a point in parameter space that is unreachable 
by any solution that starts from the origin. 

First we determine the asymptotically correct starting 
point by considering solutions corresponding to regular 
metrics at the origin. Thus we require ^(0), n(0), p(0) 
and v{0) be finite and p'(0) — 0. Moreover, the radial 
coordinate can be scaled so that A{Qi) = 1 as usual. Ex- 
panding around the origin and solving the field equations 
gives, up to leading order: 
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r\ (8) 



where the subscript refers to the values at the origin. 
Pressure and density are constant at the level of the ap- 
proximation. Regularity of the solution at origin leaves 
three free parameters vq, po and po, which completely de- 
termine the solution and hence the stellar structure. For 



a given equation of state p = p{p) the number of parame- 
ters further reduces to two. Using these equations allows 
one to start the numerical integration a small distance 
from the origin while preserving the correct asymptotic 
behaviour. 

When one relaxes the requirement that the metric and 
f (r) are to be finite at the origin, it is possible to find 
several mathematically plausible solutions at the vicinity 
of the origin with finite density p(0) but having singular 
behaviour of metric. As we will see, such solutions arise 
easily whenever one tries to solve the field equations in- 
wards, from the stellar boundary to the center. 



C. Parametrisation of density and pressure 

In order to consider more realistic matter distribution 
than predetermined toy models, we consider here poly- 
tropic stars i.e. stars with an equation of state p = np^ . 
Here k and 7 = 1 + 1/n are constants and n is often 
referred to as the polytropic index. Such equations of 
state are useful in studying white dwarfs, neutron stars 
and can also be used as a simple model of main sequence 
stars such as the Sun [3^. For a polytropic equation 
of state, one can straightforwardly solve the continuity 
equation, Eq. ([5]): 



p{r) 



= k1/(1-7) 



:^N(l-7)/(27) 

B ' 



1 



1/(7-1) 



(9) 



Requiring, that p vanishes at the stellar surface r = rn, 
where B — Birn), sets 7 > 1. (It can be shown that 
requiring finite radius constrains 7 > 6/5 in the Newto- 
nian Lane-Emden -model [39|LJ Similarly, since p' is also 
vanishing at the boundary [36|, we must further require 
that 7 < 2. 

For numerical work, it is advantageous to use 
scaled variables, p — pq6^^'^^^^\ r = ax, a = 

^ (see eg. ^9|). In these vari- 



(7- 



v/«7/(4vrG(7-l))p^' 
ables, the fundamental equation of a Newtonian star re- 
duces to the Lane-Emden -equation. In the case of GR 
or modified gravity, this is not the case, but the same 
change of variables is still useful. Using 6 instead of p is 
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advantageous also due to the fact that unhke p', 9' does 
not vanish at the boundary of the star, making identifi- 
cation of the star's surface easier. 



III. NUMERICAL RESULTS 

Our next task is to compute numerically stellar pro- 
files for certain polytropic cases for the f[R) = R — ^"^ /R 
model. We have done this for a number of polytropic 
equations of state, both starting from the center of the 
star and from the boundary. The computed profiles 9 
as well as metric components, A and n, can then be 
compared to corresponding functions of Newtonian poly- 
tropes determined by the Lane-Emden -equation. Note 
that we require that stellar solutions have finite radii, 
unlike e.g. toy models where the profile is approximated 
by an exponential function. 

We consider the Sun, as a representative of main se- 
quence stars, white dwarfs as well as neutron stars. 
The polytropic model of the Sun, with 70 — 1.2985, 
po = 1-53 X lO^kgm"^ and po = 3.00 x IQiSNm""* gives 
a fair approximation to the Stellar Standard Model [4^ . 
For the relativistic white dwarfs, we use 7 = 4/3, k — 
1.4 X 10"'^(kgm"^)-i/3 and for nonrelativistic neutron 
stars 7 = 5/3, k = 3.5 x 10-"(kg m-3)-2/3[3|^ 



A. Solution starting from the center outwards 

Using Eqs ([8]) as a starting point, we can now integrate 
numerically the field equations along with the structural 
equation, Eq. ([7]), for a given stellar model. Integration 
is stopped at the stars surface i.e. when 9 = Q. Fixing 
the central density, poj fixes also po via the polytropic 
equation of state but remains as a free parameter. In 
Fig. [l]we show the density profile for the Sun for a range 
of values of vq. We see how changing vo scales the profile 
so that a larger value of Vq leads to a star with a larger 
radius and vice versa. The corresponding Lane-Emden 
-solution i.e. the polytropic model of the Sun is equal to 
choosing uq = with very high precision. In Fig. [2] we 
show the evolution of w as a function of the scaled radial 
distance x. Two properties are notable: the evolution 
is very small i.e. the value changes very little over the 
radius of the star and the value is monotonically decreas- 
ing. The latter property is important when we consider 
fitting the star to a SdS-spacetime {v = 0) outside the 
star. 

Although the density profile in general closely resem- 
bles the Newtonian one, the behaviour of the metric is 
completely changed. Fitting the general metric to the 
PPN SdS-solution, 
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FIG. 1: Density profiles of the polytropic model of the Sun 
for diflterent values of vq. 




FIG. 2: Evolution of v{r) 
0, r = ax. 



F{r) - 4/3 for the Sun, vq 



we can solve the PPM-parameter "fppN at distance r as 



JPPN 



B(r) 
rB'(r) 



1 



Mr) 
rA'{r) 



(11) 



(10) 



Note that here we have neglected the H^r^ -term as the 
cosmological parameter is extremely tiny compared to 
relevant solar system scales. 

Numerical work shows, that at the stellar boundary 
rp jppN tends to be near value jppN = 0.5, with small 
variations depending on the value vq. In Fig. [3] we show 
the evolution of jppN outside the Sun for the vq — 
case but other choices of uq give essentially identical re- 
sults. The corresponding physical distance is equal to 
r w X X 8 X 10* m. From the figure we see that the space- 
time outside the Sun in the CDTT-model is in gross vio- 
lation of the experiments, remaining very close to 1/2 far 
outside the Solar System. The behaviour at very large 
X is explained by the fact that we ignore the cosmolog- 
ical term in the calculation of ^ppN, Eq. pT|) . The 
effect of the cosmological term should be included when 
2GM/r ^ H^r"^ ^ M^r-^, i-e. x ^ 10^\ in good agree- 
ment with the numerical calculation. 
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FIG. 3; Evolution of the PPN-parameter 'yppN outside the 
Sun. 



where = /i^/(4-\/3). At the boundary we require that 
n{rji) = B'^^^{rR)/Bext{rR), A{rii) = Aextirn) as well as 
vii^B.) = v'{rji) = 9{r}i) — 0, corresponding to the SdS- 
metric jscj . In practice we first fix central density p(0), 
and then set the gravitating mass M and stellar radius 
rR, by solving equations from the center outwards for a 
given vq. The computed mass and radius are then used 
as parameters in the outside metric, which then fixes the 
boundary condition for n and A. In Fig. IH^a), we show 
the stellar profiles corresponding to different choices of 
Vq i.e. corresponding to different values of M and r/j. 
From the figure we see that for vq = 0, or the Sun, the 
density profile is unchanged i.e. fixing the outside solu- 
tion to the SdS-metric leads to a physically acceptable 
density profile. When vq deviates from zero, the solution 
diverges near the origin. In the same figure we also show 



Note that the scale IppN at which we expect the value 
of 7 to approach unity, i.e. the range of the effec- 
tive scalar in corresponding scalar-tensor theory (see e.g. 
[13, m, US and reference therein), is IppN ^ ~ 
10^^ which far beyond the scale at which the cosmolog- 
ical term becomes effective. In the context of this par- 
ticular f{R) model, if one wishes to have jppN — > 1 on 
scales where we can ignore the cosmological term, i.e. 
2GM/r ^ /i^r^ when r ^ sets /i > 1/km for the 
Sun, M « 2 X lO'^" kg. This shows how in this model we 
cannot choose the different relevant scales to have phys- 
ically meaningful values. 

We have performed corresponding calculations also 
for relativistic white dwarfs and nonrelativistic neutron 
stars. Although stellar numerical values and scales are 
clearly different, qualitative conclusions remain: density 
profiles are close to Newtonian ones whereas metric pro- 
files are completely different and jppn ~ 0.5. We have 
also considered configurations with pressure comparable 
to density in which case we find that "fppN can deviate 
significantly from 1/2. 



B. Solutions with the SdS metric as the external 
solution 

We have seen that starting with a regular solution 
at the origin leads to unacceptable space-time out- 
side the star, i n g ood agreement with previous work 
[27I m, m, [33, H^- This result suggests that starting 
with a physically acceptable outside solution, the SdS- 
metric, will lead to an irregular solution at the center of 
the star. 

In order to determine the relevance of the SdS-metric 
for the particular f{R) model in question, we solve Eqs 
([5]) starting from the stellar boundary at given radius 
rp. We set the boundary conditions by requiring that 
external metric is the SdS-metric 
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(a) Density profiles. 




i?e.*(r) = A^xtir)-' = 1 - ^ - H^r^ 
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(b) Profiles of v{r), r = ax. 
FIG. 4: Stellar configurations with external SdS-metric. 



the corresponding evolution of v, Fig. |l{b). In all of the 
cases the conclusion is the same: v diverges at the ori- 
gin. This is as expected since we have already seen that 
any solution with a regular metric and density profile at 
the origin leads to a experimentally unacceptable outside 
solution. 

Requiring the SdS-solution as the outside metric leads 
to a divergent v and hence the scalar curvature diverges 
at the origin. Similarly also the metric components A and 
n, are irregular at the origin. In summary, if one relaxes 
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the requirement of regularity of the metric at the origin, 
one can have the SdS-metric as the outside solution and 
reproduce the density profile. 

It is worth noting, that if the star does not follow poly- 
tropic equation of state, especially in the core of the star, 
non-singular solution with external SdS-metrics can be 
found when vq = 0. In particular the regularity of the 
scalar curvature is seen directly from Eq. ^ if the equa- 
tion of state near the center of the star is relativistic 
p = p/3. We have confirmed this phenomenon by using 
an equation of state, which is polytropic in the outer re- 
gion and relativistic in the core. This resembles the case 
of massive, relativistic neutron stars fsS. 



IV. COMPARISON WITH GENERAL 
RELATIVITY 
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FIG. 5: Comparison of the size of different terms of the field 
equations for the Sun. 



We have seen that the case vq = reproduces the result 
from using the Lane-Emden -equation very well. We can 
understand this behaviour analytically by considering the 
field equations in GR and in the f(R) theory. 

In GR, the field equations can be written in the form 
R"^ = 8ttG{TI^ - where T = Ti^. In the limit 

of negligible pressure, the 00-component reads as Rq ~ 
AirGp. The scalar curvature, R, follows the density i.e. 
R^ p. 

In the CDTT- model, the situation is different. Now, as 
numerical calculations show, inside and outside the star, 
f{R) ~ i? ~ /i^ and F ^ 0(1). From the field equations 
one would then expect that for the 00-component, the 
Fi?n and OF terms to be dominant i.e. 
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SnGp or if one considers the contracted equation, Eq. 
da]), SDF « SnGp. In Fig. [5] we plot the relevant terms 
for the Sun (here Vq = 0, but the situation is unchanged 
if Vq is varied as long as 3> /). From the figure we see 
that inside the star the approximations hold and hence, 
we can write FRq « IGirGp/S. From the numerical work 
we furthermore know, that the value of u or F changes 
very little inside the star and hence we can approximate 
F « const. = Fq = 4/3 + vq. We have then 
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(13) 



Comparing to the GR expression, we see that there ex- 
ists a simple scaling between the two expressions. In the 
Lane-Emden -equation this scaling signals its presence in 
the scaling of the radial coordinate, r = ax cx x/^/G, if 
po is fixed i.e. we expect that 
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In Fig. [6] we show the Sun's profile with vq = 1 (solid 
red line) along with the Lane-Emden -solution 9le{x) 
(dotted green line) and a scaled Lane-Emden -solution, 
Ole{x/ a/I -|- (3/4)tio). The line at the bottom is the dif- 
ference between the scaled Lane-Emden -solution and the 



numerical solution for vq = 1 magnified by a factor of 
10^. As we can see, the scaled solution reproduces the 
numerical solution very well. 
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FIG. 6: Comparison of exact /(_R)-model, Lane Emden and 
scaled Lane-Emden -model density profiles. 



One can hypothesize that a similar scaling would exist 
in a general f{R) theory, f{R) ~ R + g{R,p). In order 
for such a theory to explain late time acceleration, the 
extra terms will in general have a new scale /i associated 
with the value of the present Hubble parameter (in prin- 
ciple one could also consider very finely tuned theories, 
where different parameters with different scales would 
conspire to produce current acceleration). The vacuum 
state of the theory is such that it possesses non-zero con- 
stant curvature, Rq ^ p? . Examples of such theories 
are /(i?) = R- p^/R and f{R) = R- p^ / R + R^/(3'^. 
In the latter case, choosing fi ^ p helps to avoid Solar 
System constraints [1^ by making the mass of effective 
scalar large. If R remains close to the cosmological value 
inside the star, like in the CDTT-model, it is natural to 
expect that / ~ i? ~ /x^ ~ Hq and F ^ 0{1) and hence 
a similar scaling property should apply. 
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V. CONCLUSIONS 

In the present paper we have analyzed the proper- 
ties of polytropic stars in a generahzed gravity model. 
In particular we have considered the f{R) model with 
f{R) = R — fJ.'^/R with the conclusion that the density 
profiles in general resemble the Newtonian Lane-Emden 
-solutions. Requiring that stellar solution is regular at 
the origin, we found that slightly varying the central 
curvature, vq, the stellar mass and radius are changed 
but preserve their functional Lane-Emden -form. How- 
ever, the metric components are drastically different from 
the Lane-Emden -case and therefore our results for the 
external metric conform to previously calculated results 
for completely pressure-less matter [s^] . In particu- 
lar the PPN parameter ^ppN outside the start is near 
IPPN = 1/2. 

If we do not require complete regularity of stellar solu- 
tion at the origin, but assume the external SdS-solution, 
we still find stellar profiles in good agreement with the 
Lane-Emden solution. Differences appear only near the 
center of the star and these deviations depend on how 
much the mass and radius differ from the corresponding 
Newtonian configuration. The interior solution for met- 
ric components and curvature for such stars are always 
singular, although stars with relativistic matter at the 
core may evade this property. 

Consequently, the f{R) = R — n'^/R model is not ex- 
perimentally suitable to describe the space time around 
the Sun. A possible way out is to relax the requirements 
set for the central boundary conditions, but a more plau- 
sible approach is to modify the functional form of f{R). 
The form of the action function f{R) should differentiate 
the cosmological Hubble scale Rq ^ Hq determined by 
2/(i?o) = Rof'{Ro) and the effective scalar mass scale 
oc l//"(i?o)- Then it may be possible to have "fppN 1 
at a distance IppN small compared to solar system dis- 
tances, possibly redeeming some of the f{R) models. 



This is, however, not possible in the f{R) = R — IJ.'^/R 
model, because both scales are controlled by a single pa- 
rameter /i, which when set to the cosmologically relevant 
value fi ~ Hq, leads to IppN ^ ^/Hq. In other words, 
the asymptotic SdS-metric is never reached. 

Seeking possible ways to save f{R) models note also, 
that as discussed in [36], a suitable choice of f{R) may 
exactly reproduce GR/Newtonian density profiles chang- 
ing only A and v. In this case we do not, however, know 
much about their behaviour neither outside the star nor 
near the center of the star and even cosmological con- 
straints are unknown. 

Possibly several more general /(i?) models are physi- 
cally acceptable, but in particular the model with f{R) = 
R — fi'^/R + R^ / 0^ may do after considerable fine-tuning 
of (3. Note however, that in this case an important dis- 
tinction compared to the CDTT-model applies. In the 
CDTT-model, Fg ~ 4/3 or vo ~ 0, inside and outside 
the star so that the field equations are effectively similar 
to the GR counteparts, Eq. (fT^. When the R^ term is 
added to the action, Fq 7^ 4/3 outside the star and hence 
if we still wish to have « 4/3 inside star in order to 
reproduce the Lane-Emden -solution, F must evolve sig- 
nificantly over the radius of the star. Otherwise, if Fq 
remains approximately constant, either the SdS-solution 
is not the correct outside solution, or the radius of the 
star will be different than in GR. This argument poten- 
tially offers a new, general constraint on f{R) theories of 
gravity, motivating further work. 
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